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Abstract 

An obstacle affecting any proposal for a topological quantum computer based 
on Ising anyons is that quasiparticle braiding can only implement a finite (non- 
universal) set of quantum operations. The computational power of this restricted 
set of operations (often called stabilizer operations) has been studied in quantum 
information theory, and it is known that no quantum-computational advantage can 
be obtained without the help of an additional non-stabilizer operation. Similarly, 
a bipartite two-qubit system based on Ising anyons cannot exhibit non-locality (in 
the sense of violating a Bell inequality) when only topologically protected stabilizer 
operations are performed. To produce correlations that cannot be described by a 
local hidden variable model again requires the use of a non-stabilizer operation. 
Using geometric techniques, we relate the sets of operations that enable universal 
quantum computing (UQC) with those that enable violation of a Bell inequality. 
Motivated by the fact that non-stabilizer operations are expected to be highly im- 
perfect, our aim is to provide a benchmark for identifying UQC-enabling operations 
that is both experimentally practical and conceptually simple. We show that any 
(noisy) single-qubit non-stabilizer operation that, together with perfect stabilizer 
operations, enables violation of the simplest two-qubit Bell inequality can also be 
used to enable UQC. This benchmarking requires finding the expectation values of 
two distinct Pauli measurements on each qubit of a bipartite system. 

1. Introduction and Definitions 

A topological quantum computer (TQC) would allow quantum information to 
be stored in quasiparticles and manipulated by quasiparticle braiding, in a way that 
is inherently robust against local perturbations (noise). Many of the most promis- 
ing candidates for experimental implementation of a useful TQC are those which 
support Ising anyons, e.g., fractional quantum Hall effect [1, 2], Majorana wires [3], 
p + ip superconductors [4] and the Kitaev honeycomb model (in the presence of a 
magnetic field) [5]. While there are large differences in the underlying physics of 
all these systems, the non-Abelian braiding statistics are such that they all possess 
the same computational power (i.e., they enable the same types of quantum gates 
and measurements when viewed as a quantum information processing device). 

The model of computation (sometimes known as the Clifford computer model 
[6]) that we assume in this work, and the one that is relevant to Ising anyon TQCs, 
is one in which (i) all single-qubit Clifford gates (discussed in Section 1.1), (it) 
measurements in the computational basis, and (Hi) the controlled-NOT gate, are 
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all implemented in a topologically protected (i.e., effectively perfect) manner. Col- 
lectively, we refer to this set of gates and measurements as stabilizer operations. 
Some proposals for TQCs based on Ising anyons use "lag) -distillation" to achieve 
the controllcd-NOT [7], while others achieve it using quasiparticle braiding [8] or 
non-demolitional measurements of the collective charge of four anyons [9]. It is 
the latter, topologically protected, implementations of the controlled-NOT that we 
have in mind here. 

A result in quantum information theory - the Gottesman-Knill theorem (see 
e.g. [10] for a discussion) - says that the set of stabilizer operations is insufficient 
for achieving universal quantum computation. If a (noiseless) single-qubit unitary 
gate from outside the Clifford group was implementable then we would immediately 
have full UQC (this is the case for TQC proposals based on Fibonacci anyons). In 
fact, for an Ising anyon TQC, attempting to implement any non-stabilizer operation 
necessitates using non-topological operations which are expected to be highly noisy. 
For specific target gates and noise models, one can calculate the threshold noise 
rate [11, 12, 13] before the power to provide UQC is lost (see e.g., Section 4.1 
for a relevant example). Rather than consider a number of different target gates 
(unitaries) and noise models, we examine an overall quantum operation £ , which we 
intend to be non-stabilizer but which may be subject to some unknown evolution, 
and give an operational benchmark on the utility of £ for UQC. Proving that 
an operation, £, enables UQC, when used alongside perfect stabilizer operations, 
reduces to the question of whether £ can be used to produce ancillas that are 
suitable for magic state distillation [14, 15] (MSD). MSD is a subroutine used in 
many proposals for fault-tolerant UQC, whereby the requisite non-Clifford gate, 
U, is implemented with the assistance of a non-stabilizer pure qubit state \ipu}- 
Crucially, many impure (noisy) copies of \ipu) can be used to create a smaller 
number of purer \ipu) using only stabilizer operations, and this process can be 
iterated a number of times if necessary. However, if the ancillas \ipu) are too 
impure then it becomes futile to attempt an approximation of U, regardless of the 
number of ancillas and techniques used. A summary of the current methods and 
limits of MSD can be found in [16, 17, 18]. 

The benchmark (for identifying £ that enable UQC) that we have just described 
will involve the violation of a Bell inequality by a two-qubit quantum state, wherein 
the underlying qubit encoding is provided for by non-abelian Ising anyons. The 
topological phase supporting Ising anyons is inherently nonlocal, but for the pur- 
poses of quantum information processing we are more interested in displaying the 
nonlocality of encoded qubits. Investigations of nonlocality in non-abelian anyons 
have been carried out in [19, 20] from a different perspective. The results here 
are motivated mainly by Ising anyon proposals for a TQC, but we note that any 
system in which magic state distillation is used to achieve fault-tolerant quantum 
computation assumes that stabilizer operations are effectively perfect (albeit at an 
encoded level within a CSS error-correcting code, for example). Consequently, the 
results presented here can be applied to any such system. 

Given the restriction on topologically implementable gates, and hence on mea- 
surement directions, one quickly sees (see e.g., Section 1.3 later) that this setup 
precludes the possibility of performing measurements that lead to violation of a 
Bell inequality. In other words, all measurement outcomes can be described by 
a local hidden variable (LHV) theory, when we restrict the experimenter to using 
only topologically protected gates. Previously, we have motivated examining the 
set of operations £ that can provide UQC. Similarly, we can examine the set of op- 
erations £ that can enable violation of a Bell inequality, when used in conjunction 
with topological operations. Here we relate the two sets of operations with each 
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other, and show that, under certain restrictions, one is strictly included within the 
other. More generally, we show that any operation that can be used to violate a 
so-called CHSH Bell inequality can also be used to provide UQC. 

Note that violation of Bell inequalities has been related to quantum computa- 
tional power in other contexts [21, 22], although these works were more focused on 
generalized entanglement and measurement-based quantum computation respec- 
tively. While an interesting topic in its own right, our results do not claim to relate 
nonlocality and UQC in a broad sense. Indeed in the Ising anyon context, if we are 
allowed to construct tripartite states (e.g., the GHZ state (|000) + 1 1 1 1 ) ) / >/2 ) then 
topologically protected operations are sufficient to exhibit tripartite nonlocality 
(e.g. Mermin's experiment [23]), whilst still not enabling UQC. 

1.1. Clifford Operations. The single qubit Clifford group has 24 distinct ele- 
ments and contains the Pauli group (of order 4) as a subgroup (strictly speaking 
we are discussing the Clifford and Pauli groups modulo their center, which amounts 
to ignoring a global phase in the matrix representation of these elements). Canon- 
ical, non-Pauli, elements of the Clifford group are 

HjC-0- s -(i?) 

and these two operations are actually sufficient to generate the whole group. The 
characteristic property of Clifford gates, C, is that they map Pauli operators to 
Pauli operators under conjugation, i.e., 

CdjC^ H> ±a k j,ke{x,y,z} 

When visualized as operations on the Bloch sphere, the Clifford group can be iden- 
tified with the 24 possible transformations that can be constructed by combining 
consecutive 90° rotations about the x-,y- and z- axes. 

1.2. Quantum Operations. To describe the unknown evolution of an Ising anyon 
qubit, whilst we attempt to implement the requisite non-stabilizer operation, we 
use the quantum operations formalism (as described in e.g. [10]), where £ is a 
superoperator that maps input density matrices to output density matrices: 

£(Pin) — Pout- 

A well known tool in quantum information, the Jamiolkowski isomorphism, tells us 
that £ is completely characterized by the output state, g$ , of a process whereby £ 
is applied to one half of a maximally entangled pair 

(2) g £ =(I®£)[|*X*|] where $ = ^Mii) 

Because gs is our preferred method for representing general quantum operations, 
we will sometimes use £ and gg interchangeably in latter sections of the paper. 

Another convenient way of expressing a quantum operation is in terms of its 
so-called Kraus operators, {Ei}, via 

£{pin) = Pout = ^2 E iPin E l (where ^ E l E i = J ) 

i i 

It takes 12 real parameters to completely characterize an arbitrary completely- 
positive trace-preserving operation. A 9-parameter subset of the set of all possible 
operations is given by those that preserve the identity (these are often called unital 
channels) i.e., those for which £(T) = I. For unital operations, a Kraus- like descrip- 
tion is possible, except the operation is now a probabilistic mixture of unitarics 
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(where each unitary Uk is applied with probability pk) 

£(Pin) = Pout = ^PiUiPinU} (where y^gj = 1). 

i i 

1.3. Nonlocal correlations with restricted operations. Here we briefly mo- 
tivate why a non-stabilizer operation is necessary to exhibit non-locality. Consider 
the maximally entangled state |<T>) = l 00 ^ 11 ^ , then its Pauli expectation values, 
defined as 

JK = Tr(|$X$|(7j<g)o-fe) j,k£{x,y,z} 
are all zero except for 

XX = -YY = ZZ = 1 

The Pauli expectation values for |<£>) can be recreated exactly by two spatially sepa- 
rated parties A and B, given that they share 3 random (unbiased) bits {r*i, r2, ^3} 6 
{0, 1}, and they both obey the following set of rules [24] 

(1) Measurement in the X direction ■<-» A and B both output — l ri 

(2) Measurement in the Y direction ■<-» A outputs — l' 1 * 2 , B outputs — l' r2+1 

(3) Measurement in the Z direction f-> A and B both output — V 3 

A little thought shows that a similar scheme would work for any bipartite entangled 
state created by stabilizer operations. That such a scheme suffices, using only shared 
randomness, indicates that no purely quantum mechanical effects are needed to 
describe such an experiment. 

2. Bell Inequalities 

We will attempt to perform an experiment that, unlike the example of Section 
1.3, is not describable by local hidden variables. The basic setup is depicted in Fig- 
ure 1, where £ is applied to one half of a maximally entangled state (|00) + |ll))/\/2 
and then Pauli measurements are performed. We can allow for either two or three 
distinct Pauli measurements to be performed on each qubit (Figure 1 depicts the 
scenario in which only two measurements are performed on each side), and these 
lead to different types of Bell inequality. 




Figure 1. A simple set-up to detect nonlocality: Here, two pos- 
sible measurement settings for the first (second) qubit are denoted 
Ai (Bj). When A; and Bj are constrained to be Pauli operators, 
then £ must be a non-stabilizer operation if any Bell inequality 
(e.g., that given in (5)) is to be violated. 
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For a given bipartite state g, and using a notation similar to that of Collins and 
Gisin [26] , we can arrange the expectation values for a given nonlocality experiment 
with Pauli measurements in a table i.e., 



(3) 



/ II XI YI ZI ^ 

IX XX YX ZX 

IY XY YY ZY 

\ IZ XZ YZ ZZ J 



where XY is the expectation value Tt(q<t x ® a y ) and so on. Since each Pauli op- 
erator has eigenvalues ±1, we can list all possible matrices of the form (3) that 
correspond to deterministic local configurations. A deterministic local configura- 
tion is one in which the local Pauli expectation values are extremal and, when 
multiplied, provide the nonlocal Pauli expectation value (e.g., where identities like 
(XI)(IY)=XY hold, with XI,IY,XY all either ±1). By letting a, &,...,/ € {0, 1} 
take on all 2 6 assignments, then the following matrix fully describes all possible 
local configurations, 



(4) 



1 (-1)" 
(_l)d 

(-l) e (-l) a+e 

V (-1)/ 



(-i) b 
(-i) b+d 

(_ l)W-e 

(_!)*+/ 



(-l) c \ 

(-iy+ d 

(_l)c+e 
(-1) C+/ / 



If an experimentalist's measured expectation values correspond exactly to any of 
the above local configurations, then these expectation values do not exhibit any 
nonlocality. Clearly, any realistic experiment will not obey such a strict condition 
and the question of exhibiting nonlocality becomes richer. The experimentalist 
must now check that the measured expectation values cannot be expressed as a 
probabilistic combination of local configurations. This naturally leads to the notion 
of convex geometry and bounded polyhedra in higher dimensions (these are usually 
called polytopes). 

Ignoring the constant term in the top- left, the 2 6 matrices in (4) can each be 
identified with a vector in R 15 . These local configuration vectors are called vertices 
and the convex hull of these vertices (i.e. the set of vectors that is expressible 
as a probabilistic combination of the vertices) describes a polytope in R 15 (this is 
analagous to how the 8 vertices (±1,±1,±1) describe a solid cube in R 3 ). The 
interior of this polytope describes all possible LHV models. Conversely, any table 
of expectation values (3) that does not lie inside this polytope exhibits genuine 
quantum nonlocality. The usual, textbook, way of identifying nonlocality is via 
Bell inequalities like 



(5) 



XX + XY + YX - YY < 2. 



In fact, (tight) Bell inequalities such as these are actually the defining inequalities 
for the bounding faces (facets) of the LHV polytope that we have just described. 
Software such as Avis' Irs [27] can be used to derive all the facets of LHV polytope, 
given the 2 6 local configuration vectors in (4) as input (and vice versa) . In order to 
unify our notation, we rearrange the inequality (5) and then rewrite in table form 



(6) 
(7) 



2 - XX - XY - YX 



2 

-1 -1 

-1 1 








o / 



YY > 



II 
IX 
IY 
IZ 



XI 

XX 
XY 
XZ 



YI 

YX 
YY 
YZ 



ZI 

ZX 
ZY 
ZZ 



> 
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where the dot product between two matrices M, N behaves like the familiar vector 
dot product i.e., M ■ N — Tr(M T N) = M i,j N i,j- The ful1 list of bounding 
inequalities for the LHV polytope is described in a concise form in Section 4. 



3. Non-stabilizer operations 

In the previous section, ascertaining whether measurement results had an LHV 
model amounted to the question of whether the results could be expressed as a 
probabilistic combination of local configurations. In the context of UQC, the rel- 
evant question is now whether £ is expressible as a probabilistic combination of 
Clifford gates; for example, an operation £ given by 

£{p) = iffpfft + * ,s p st + $_ p 

is clearly of no use in the quest to achieve UQC. The general prerequisite for £ to 
be useful is 

24 / 24 \ 

(8) £ (p) ? ViGipC\ where £ Pi = 1 

i=l \ i=l / 

If we use the Jamiolkowski isomorphism (as defined in (2)) to represent operations, 
then a completely equivalent condition to (8) is given by 

24 / 24 

(9) Q £ ^Ypi\J Ci ){J Ci \ where 

i=l \ i=l / 

To test whether an unknown operation satisfies the above requirement, one is nat- 
urally lead to the concept of convex polytopes once more. The 24 Clifford gates, 
when represented as \Jci)(Jci\ and decomposed in the Pauli basis, form the 24 
vertices of what is called the Clifford polytope. In fact, the Clifford operations only 
span a 9-dimensional subspace of R 15 because the six coefficients IX,IY,IZ,XI,YI 
and ZI are all identically zero when £(T) = I. The Clifford gates given in (1.1), for 
example, correspond to the following two vertices 



^Tpi = 1 and \J Ci ) = (I ® d) 



( 



\ 



v s 



( 



\ 



1 / 



The facets of the Clifford polytope [11, 13, 28] are bounding inequalities that 
partition the set of all valid qs into those that represent £ satisfying (8) and those 
that do not. In other words, any operation £ that satisfies (8) violates (at least) 
one of these facet inequalities. These inequalities, which we label /, can be cast in 
matrix form, in a completely equivalent way to what was done in (7), and it turns 
out there are two distinct classes of facet [11] - a representative example of each of 
which wc provide below. 



(10) I a = {all facets of type I a or /J}, 



(11) I p = {all facets of type Jg}, 



The total number of facets, required to completely characterize the boundary of the 
Clifford polytope, is 120, which can be broken down into |X a | + |Tys| = (24 + 24) + 72. 
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-1 
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We are primarily interested in facets from Ip because it was shown in [13] that an 
operation violating such a facet can always be used to create an ancilla p that is 
useful for magic state distillation. The circuit to create p is given in Figure 2, where 
the choice of which two-qubit Pauli measurement II to perform is determined by 
the particular facet that is violated by £. 



— H 




n 






e 






5- 





-|0) 



Figure 2. A circuit to help achieve universal quantum compu- 
tation via magic state distillation: Every element of this circuit 
except £ is implementable using stabilizer operations. When £ 
exhibits nonlocality in the setup of Figure 1, then £ used in the 
above circuit produces ancillas p that are useful for magic state 
distillation subroutine (MSD circuit not depicted). The block con- 
taining II stands for a two-qubit Pauli measurement (e.g. parity 
measurement) wherein we postselect on the desired outcome. 



4. Results 



A full facet description of the LHV polytope as described in Section 2 comprises 
684 distinct facets, and this facet description can partitioned into 3 different classes 
i.e, 



(12) l tnv = {all facets of type Itnv}, Itriv = 



(13) I2222 = {all facets of type 12222}, ^2222 = 



(14) I3322 = {all facets of type J3322}, ^3322 



V 

4 

-1 
-1 




1 

-1 



2 








-1 
-1 











-1 
1 





-1 
1 
1 

-1 



-1 
1 
1 
1 






\ 















/ 





\ 















/ 





\ 


-1 




1 







/ 



The number of distinct facets in each class is \T tr iv\ = 36, IX2222I = 72 and 
IX3322 1 = 576. Inequalities of the type I tr iv are considered trivial, since the restric- 
tions they impose are satisfied by any bipartite quantum state g. The representative 
inequality Itriv given in (12) amounts to Tr [(I — a x ) ® (I — <j x )q] > 0, for example. 
Inequalities of the type ^2222 are usually known as Clauser-Horne-Shimony-Holt 
(CHSH) inequalities [29] and an example of one such inequality was already dis- 
cussed in the context of Equation (7). These -Z2222 inequalities appear, in a different 
form closer to that of (5), in most textbook accounts of nonlocality (see e.g. [10]). 
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Bell inequalities of the type -Z3322 are less well known but have a long history, ap- 
pearing as early as 1981 [30] and rediscovered more recently in [31, 32, 26]. Finally 
we note the reasoning behind the notation I cc dd', the numbers cedd in the subscript 
refer to the number of measurement settings, c, available to each party, wherein 
each measurement has a number, d, of possible outcomes (e.g., inequalities of the 
type -Z3322 can only arise when there are 3 possible measurement settings on each 
side of the bipartite state). 

The -Z2222 inequalities require only two Pauli measurements on each side whereas 
^3322 inequalities require three on each side. It turns out that the ^2222 inequalities 
are sufficient for our purposes. In fact, ^3322 inequalities are irrelevant (redundant) 
for a very large subset of all possible quantum operations £. If one makes the 
assumption that £(T) = I, which is the case for many of the most important noise 
models, then local expectation values for g are identically zero, and coefficients 
of IX,IY,IZ,XI,YI and ZI in 13322 facets can be ignored. With this modification, 
the -Z3322 inequalities are expressible as linear combinations of -Z2222 (CHSH) in- 
equalities, and hence the 13322 inequalities are irrelevant. One can check that all 
(purely nonlocal) ^3322 inequalities decompose into four J2222 inequalities, as in the 
following example 



(15) 



\ 



-1 
1 






( 2 
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i 2 
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-1 
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1 












1 
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+ 
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/ 
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-1 


1 





) 



We are now ready to state the main result: 

Theorem 1. If 'there is an operations violating a facet from I2222 (13) (i.e., violat- 
ing a CHSH inequality), then £ also violates a facet (of the Clifford polytope) from 
Tp (11). Using results from [13], such £ enable UQC (via magic state distillation) 
when supplemented with topologically protected operations. 

Mathematically, we have the following statement 

(16) V physically valid £ , V ^2222 € ^2222, 3 Ip 6 Ip such that 

violates Ip at least as much as it violates ^2222 

Without loss of generality, we assume that the canonical representative, ^2222 (13), 
of X2222 is violated, and show that Ip (11) homTp is necessarily also violated. This 
is a consequence of the following proposition, which is easily proved 



(17) 



1000 
0-1-10 
0-110 
0001 



II XI YI ZI 
IX XX YX ZX 
IY XY YY ZY 
IZ XZ YZ ZZ 



< 



2 
0-1-10 



-1 




Rearranging the above inequality gives 



(18) 



V 



-1 







\ 



II 


XI 


YI 


ZI \ 


IX 


XX 


YX 


ZX 


IY 


XY 


YY 


ZY 


IZ 


XZ 


YZ 


ZZ / 



II XI YI ZI 
IX XX YX ZX 
IY XY YY ZY 
IZ XZ YZ ZZ 



< 



which simply stands for 
(19) 



-1 + ZZ < 0. 



This is always satisfied since Pauli expectation values are in the range [— 1 , 1] . 
create an ancilla for magic state distillation (and consequently enable UQC), 



To 
the 
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best measurement II to use in the circuit of Figure 2 is the parity measurement 
postselcted on even parity i.e., II = \ (I + a z a z ). 

Note that = IZ2222I = 72. Examining the matrix representation of these 
sets of facets, one sees that each Ip s Ip can be paired one-to-one with each 
^2222 G I2222, by matching facets whose non-identity coefficients differ in only one 
position (e.g., the coefficient of ZZ in our canonical ^2222 (13) is 0, whereas for our 
canonical Ip (11) the coefficient is +1). This provides a recipe for identifying the 
relevant Ip such that (16) holds, for a given J2222 (the CHSH inequality that we 
presume has been violated). Any such pair {12222,^/3} will obey an inequality of 
the form (17). The general expression, analagous to (19), for different facets from 
I2222 and suitably chosen facets from Ip is 

(20) -1 ± Trfo ® a kQe ) < j,ke{x,y,z} 

and this clearly always holds. In order to use £ to obtain ancillas useful for magic 
state distillation, a suitable measurement II in Figure 2 includes one of the form 
n= ±(I±a j( Tk) [13]. 

4.1. Example: Phase Gate subject to Dephasing. The phase gate, U z (8), 
defined as 

is an element of the Clifford group when 8 is a multiple of ir/2. Whilst techni- 
cally any non-Clifford angle 8 would suffice to enable UQC, one typically seeks to 
implement the so called "pi-over-eight" gate £4(f ) because it possesses additional 
desirable properties [33]. For the case of Ising anyons, it has been proposed [34, 35] 
to perform a phase gate using a sack-like geometry, wherein the sack contains a 
single qubit, and an anyonic edge current can either follow the exterior boundary 
of the sack or tunnel across the constriction of the sack. In [34] it was shown that, 
subject to certain assumptions, the overall evolution of the qubit state can be de- 
scribed by a superoperator £(p(to)) = pit) where the matrix elements of the density 
operator change like 

u\ = f Poo(to) e-^e- tf> p Q1 (t ) \ 

or, equivalently in the Kraus operator form, 

(21) p(t) = E p(to)El + E lP (t )E\ where 

As a side note; in quantum information theory, the description of this process in 
terms of Kraus operators would probably be rewritten 

(22) £(p) = EqpEI + E lP E\ with 

E = v/T - p U z {8), E 1 =^pa z U z (6), where p = 1 - e 

which is completely equivalent. The natural interpretation is that with probability 
1 — p the desired U z (8) gate is performed, and with probability p an undesired 
rotation cr z U z (8) is performed. This undesired rotation is actually the worst possible 
noise that could be inflicted [12] (i.e., it requires the smallest noise rate, p, to make 
the overall £ expressible as a probabilistic combination of Clifford gates). 
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The operation, £, of (21), when expressed as a Jamiolkowski state, gg 
decomposed in the Pauli basis, takes the form 



and 



(23) 



V 



II 

IX 
IY 
IZ 



XI 
XX 
XY 
XZ 



YI 

YX 
YY 
YZ 



ZI 

zx 

ZY 

zz 






V o 





L Cos[0] e" 

" Sin [6»] -e 






t Sin [9] 

£ Cos [0] 







1 / 



and so the Bell inequality given in (7) reads 

XX + XY + YX - YY < 2 =>■ {[cos 9] + [sin.0] + [sin0] - [- cos0]) < 2 

For the optimal choice of angle, 9 = this inequality is violated for the parameter 
range < s < Vm2. In terms of the simplified Kraus operators given in (22) 
the range of allowable noise rates p, while still giving CHSH violation is given by 
< V < 14%. 

Another simple calculation shows that the facet Ip is also violated (hence UQC 
is possible) for the same noise rates < p < 14%. One can actually see this straight 
away by noting that (23) implies that — 1+ZZ= 0, which means (with reference to 
(19)) that the threshold noise rate is the same for violation of Ip and 72222- This will 
also be true for any noise model that is expressible as a probabilistic combination 
of rotations about the z-axis (and similarly for operations comprised of rotations 
about x and y axes too). 

The noise rates for which Bell inequality violation is possible, and the noise rates 
for which UQC is possible, generally do not coincide. For example, a U z (j) gate 
under depolarizing noise, modeled as 

£{p) = {l-p)U z {l)pU z {iy +p\ 

enables UQC, but does not violate any Bell inequality, for the parameter range 
0.29 < p < 0.45. In the current context, a non-Clifford operation £ is necessary if 
one hopes to observe nonlocality (as shown in Section 1.3), but not always sufficient. 
However, if we expand the problem to allow multiple (possibly simultaneous) uses of 
£, in a more complex circuit than that of Figure 1, then it is an interesting question 
whether £ being non-Clifford could be both a necessary and sufficient condition for 
the violation of a bipartite Bell inequality. 



5. Conclusion 

By using techniques from convex geometry, we can show that detection of non- 
locality in the encoded qubits of an Ising anyon TQC serves as a benchmark for 
universal quantum computation. The strength of this approach is that it does 
not consider specific noise models, but rather considers general quantum opera- 
tions. Given an unknown operation that has been observed to enable violation 
of a CHSH inequality (using two distinct Pauli measurements on each half of an 
entangled state), we can prescribe a way of using this operation to manufacture 
ancillas that are suitable for magic state distillation. When used in conjunction 
with topologically protected stabilizer operations, ancillas of this kind allow for 
universal quantum computation. 
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